In-event background and signal reconstruction for 
two-photon invariant-mass analyses 



O 
O 

(N 

C 



Nick van Eijndhoven and Wouter Wetzels 

Department of Subatomic Physics, Utrecht University /NIKHEF 
P.O. Box 80.000, NL-3508 TA Utrecht, The Netherlands 



Abstract 



A method is presented for the reconstruction of both the background and signal in invariant-mass 
analyses for two-photon decays. The procedure does not make use of event mixing techniques and 
as such is based exclusively on an event-by-event analysis. Consequently topological correlations of 
the event (e.g. jet structures) are automatically taken into account. By means of the decay process 
7t° — > 77 it will be demonstrated how the procedure allows for determination of the 7r° yield from the 
observed decay photons. 
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Q_i, 1 Introduction 

Oh! 

In particle physics experiments some particles have such a short lifetime that direct de- 
tection is impossible. In these cases the decay products have to be used to reconstruct the 
parent-particle distributions. For example, an invariant-mass analysis allows to determine the 
^ ■ yield and also the transverse-momentum distribution of the parent particles [1]. 

In invariant-mass analyses it is essential to separate the invariant-mass peak from the 
combinatorial background. Using an event mixing procedure an approximation for this com- 
binatorial background can be obtained. The invariant-mass peak is then reconstructed from 
the total invariant-mass distribution by subtracting this background. In order to construct 
the background by means of such an event mixing procedure, particles originating from dif- 
ferent events have to be combined in pairs. Differences in the characteristics of the events 
used in the mixing procedure will cause the constructed background to deviate from the 
actual combinatorial background. The degree of deviation from the actual background de- 
pends on the position correlations of the decay particles. For instance, in the determination 
of the transverse-momentum distribution of the parent particles, this background has to be 
constructed for various transverse-momentum values. The influence of jet structures might 
be restricted to a certain transverse-momentum range, which implies that the correlations of 
the decay particles are related to the transverse momentum of the parent particles. In this 
way a systematical effect is introduced which affects (the slope of) the transverse-momentum 
distribution. 

From the above it is seen that it would be preferable to use a method that does not 
require combination of particles from different events. Such a method, which takes topological 
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correlations of the particles within the event into account, is presented here. It is applied to 
the analysis of the decay process tt° — > 77 under circumstances matching those of Pb+Pb 
collisions at the CERN-SPS. The datasets used for evaluation of our analysis method were 
obtained by means of a computer simulation that models the 7r° spectra as observed in heavy- 
ion collisions at SPS energies [2]. 

A realistic energy resolution has been implemented to mimic the photon detection in an actual 
calorimeter system like the one of the WA98 experiment [2]. The details of this simulation 
are given in section 6. 

2 The invariant mass distribution 

Consider a heavy-ion collision in which a certain amount of 7r°s is produced. More than 
98 percent of these 7r°s will decay into a pair of photons. Because of their short lifetime 
(cr = 25.1 nm), the 7r°s can be regarded to decay immediately after their creation at the 
interaction point. Under realistic experimental conditions only a fraction of the decay photons 
will be detected. The energy and position of a photon i are recorded in the detection system 
and will be denoted by Ei and afj, respectively. Consequently, the four-momentum of each 
detected photon is determined. 

The invariant mass Mi nv of a pair of photons with four-momenta p\ and p 2 is computed in 
the following way : 

ML = (P1+P2) 2 = 2(E 1 E 2 - pi -p 2 ) = 2E 1 E 2 (1 - cos(0i 2 )) , (1) 

where $12 is the angle between the momenta of the two photons. 

All possible combinations of two photons within an event will provide a distribution con- 
sisting of the 7T° invariant-mass peak and a combinatorial background. The latter being due 
to combinations of two photons which did not both originate from the same 7r° decay. The 
distribution obtained in this way will be called the total invariant-mass distribution T. 
The combinations of energies and positions that constitute this T distribution are indicated 
below 

E\X\ ® E 2 x 2 (photon 1 combined with photon 2), 
E1X1 ® E 3 x 3 , 

E N _ix jv„i ® E N x N , 



in which a combination E^Xi <8> EjXj yields an invariant mass Mjj = ^J2EiEj(l — cos(i?jj)). 
An example of such a total invariant-mass distribution, resulting for events containing 50 
each, is shown in fig. 1. Here the tt° peak at about 135 MeV is seen positioned on top of 
a large combinatorial background. Since the integral of the peak corresponds to the number 
of 7r°s for which both decay photons were detected, being the basic ingredient in the studies 
of 7T° spectra, it is essential to separate the tt° peak from the background. 
This cannot be done with sufficient accuracy in a straightforward way, since many combina- 
tions of photons that contribute to the background have an invariant mass close to the ir° 
mass. In addition, the width of the reconstructed ir° peak is relatively large, due to the limited 
energy and position resolution of the detection system, and the shape of the combinatorial 
background is unknown. 
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If an accurate prediction of the (shape of the) background distribution could be obtained, the 
7T° peak could be retrieved from the total invariant-mass distribution T simply by subtracting 
this background. 
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Fig. 1. The total invariant-mass distribution T for the case of 50 7r°s/event. 

3 Background reconstruction via event mixing 

A distribution that accurately represents the combinatorial background can be obtained 
by means of an event mixing method. For each entry in the "event mixing" background 
distribution, two sets of photon momenta (obtained from two different events) are used. One 
photon momentum is taken from the first of these two sets and another is taken from the 
second. These two are then combined to provide an entry in the invariant-mass distribution. 
The combinations used in the case of mixing events A and B, of which the numbers of detected 
photons are Aa and Ab respectively, are given by : 

£l[Apl[A] ® £l[Bpl[B]> 
£l[Apl[A] ® E 2 [B]X2[B], 

E Na [A] X Na [A] ® E Nb [ B ]X Nb [b] • 

The invariant-mass distribution obtained in this way is expected to describe the actual 
combinatorial-background distribution rather accurately. No combinations of two photons 
originating from the same 7r° decay occur. Therefore, the 7r° peak will be absent. 
An imperfection of this method is the fact that the two events used for mixing will have differ- 
ent position and energy distributions. In addition, for heavy-ion collisions also the centrality 
of the two events is different. The latter may be partly overcome by grouping the collisions 
into centrality classes and combining only events that belong to the same class. Still, possible 
correlations in photon positions within the original event will surely be affected when a second 
event is used in the mixing procedure. 

Since in general the 7r° peak is relatively small compared to the background, the error in the 
derived number of parent particles due to these effects might become substantial, especially 
in the case of constructing transverse-momentum spectra. 
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4 Background and signal reconstruction via position swapping 

Below we present an alternative method to reconstruct both the combinatorial background 
and the invariant-mass peak signal of the parent particles. The procedure is carried out within 
one and the same event and as such does not require mixing of data from different events. 
In order to perform this analysis a new distribution has to be constructed, to be called the 
"position swapped" distribution S. For each entry in this distribution two photons are taken 
from a single event. However, before computing the invariant mass of this combination, one 
of the two photon positions is replaced by the position of a randomly chosen third photon 
from the same event. The combinations which constitute this <S distribution are given by 
EiXi <8> EjX a , where i and j indicate photon i and j as before and a denotes the randomly 
chosen third photon (a ^ i ^ j). 

Unlike in the case of event mixing, only one event is used to compute each entry. An im- 
portant advantage is that when the photon positions are swapped the global characteristics 
of the event are not modified, thus preserving possible topological correlations. 
In the case of azimuthal symmetry and a significant energy-position correlation for the pho- 
tons, it is recommended to divide the photon sample into categories characterised by the 
photon polar angle. In this case the energy-position correlation can be maintained in the 
position-swapping process, by only combining photons belonging to the same category. 

To enhance the statistical significance, the number of entries can be increased by allow- 
ing a to take multiple values. Furthermore, it is recommended to replace both Xi and Xj 
subsequently as described before to make sure that all of the photon positions occur with 
comparable statistical weights. In order to extract the signal, the S distribution has to be 
normalised properly, as outlined hereafter. 

An example of the resulting S distribution corresponding to the event sample reflected in 
fig. 1 is shown in fig. 2. 
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Fig. 2. The background distribution S for the case of 50 7r°s/event. 

It is instructive to regard the S distribution as being composed mainly of a combinatorial 
background with in addition small "matching energy" (£) and "matching position" (V) com- 
ponents. The meaning of these £ and V components will be explained hereafter. 
Consider the case that photon k and photon I are created by the decay of a certain tt°. 
As a consequence, the combination E^x^ <8) E[Xi will then contribute to the 7r° peak in the 
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total invariant-mass distribution T. In the position-swapping procedure, the combination 
E k x k <g) E\x m will be encountered in the case that the randomly chosen photon is photon m 
(m 7^ k 7^ I). However, in this particular combination the energy E\ happens to be the energy 
of the photon which actually originated together with photon k from the decay of one and 
the same 7r°. The resulting entry in the invariant-mass distribution thus contains two photons 
with what we call "matching energies" and consequently contribute to the £ component. 
In a similar way, the combination E k x k <8> E m X[ contributes to the V component in which 
both photon positions correspond to those of the two photons originating from the 7T° decay. 
Here the randomly chosen photon happens to be photon I, whereas the initial photon pair 
consisted of photons k and m. 

The £ and V components of the S distribution are intrinsically different from the combinato- 
rial background, as outlined below. The remaining combinations entering the S distribution 
yield a distribution in very good correspondance with the actual combinatorial background. 

Consider the case that for some event N p h ot photons are detected and that the ir° invariant- 
mass peak for this event contains N peak entries resulting from actual 7r° decays. The combi- 
natorial background in the total distribution T then consist of [\N p h ot {N p h ot — 1) — N peak ] 
entries. The position-swapped distribution S for this particular event has an £ component 
with exactly N peak entries. The V component will, statistically, have N peak entries as well. 
The remaining [^N p h ot (N p h ot — 1) — 2N pea k] entries in the S distribution contribute to the 
combinatorial background. 

It is crucial to accurately determine the £ and V contributions in order to determine the 7T° 
peak contents. First, the procedure to construct the £ distribution will be described. It will be 
explained hereafter that in this process an additional distribution IA is needed, consisting of 
entries of the form \J ]Zcole^ - Here 6a and 9b represent random angles between two photon 
momenta from the same event. This distribution is normalised to have an integral equal to 
unity in order to ease a proper scaling to the actual peak contents afterwards. 

The invariant mass Mg of the combination E k x k (£> Eix m contributing to the £ distribution, 
in which m is the randomly selected photon, is given by 



M £ = y/2E k Ei(l-cose km ) = \j2E k E^\ - co S 9 kl )J {1 Z^^Y = M W T[ 



1 - COS 9 kr . 



(1 - cos6> fc /) 

(2) 



As indicated by this formula, the £ component of the S distribution is obtained by combina- 
tion of the 7T° peak (W) with the U distribution in the following way : 

£{E) = du U{u)W{E w )5{E - uE w ) dE w . (3) 

JO JO 

In a similar way we can construct a distribution that represents the V component of the 
S distribution. For this a distribution V is needed that is constructed from ratios of random 

energies, i.e. \ as indicated hereafter. 

V 1 

The invariant mass M-p of the combination E k x k <8> E m x\ is : 



M V = ^2E k E m {\ - cos9 kl ) = ^2£*£,(1 ~ cos^/fr = M uJ^r- ( 4 ) 
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Consequently, the V component of the S distribution is obtained by combination of the V 
peak (W) with the V distribution in the following way : 



V{E)= dv V(v)W{E w )8(E-vE w )dE v 
Jo Jo 



(5) 



We have now determined all the necessary components to construct the background of the 
two-photon invariant-mass spectrum. This will allow us to extract the number of entries in 
the 7T° peak due to genuine 7r° decays and thus the actual number of detected 7r°s within a 
certain event sample. 

The distribution that results when the <S distribution is subtracted from the T distribution 
will be called the "difference" distribution V. It consists of the N pea k entries in the 7r° peak, 
Npeak entries for the combinatorial background, minus the £ and V components of the S 
distribution. An overview of all these contributions is provided in table 1. 
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Table 1. Composition of the various distributions in terms of the number of entries. 

An example of this T> distribution for the case of 50 7r°s/event is shown in fig. 3. It was 
obtained by subtracting the S distribution of fig. 2 from the T distribution of fig. 1. 
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Fig. 3. The V distribution for the case of 50 7r°s/event. 

We now represent the actual 7r° peak by a gaussian of which the mean, standard deviation 
and the number of entries are regarded as free parameters. The £ and V components corre- 
sponding to each of the trial gaussian distributions are obtained in the way described above. 
In this procedure the position-swapped distribution S provides an accurate description of the 
combinatorial background. The four distributions T, S, £ and V are then combined to give 
a prediction for the T> distribution, denoted by T> W ed- 
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In order to obtain the best matching of the calculated (T> pre( i) with the actual T> distribution, 
a x 2 value per bin is computed for each of the trial gaussian distributions. This is performed 
in the 2a interval around the 7r° mass. The reason that we restrict ourselves to the 2a region 
for the fitting procedure is that in the tails the actual ir° peak is not accurately described 
by a gaussian distribution, which would have a degrading effect in the fitting procedure. 
The gaussian distribution that produces the T> pre d distribution that fits best to the actual T> 
distribution represents the actual ir° peak in the interval mentioned above. Consequently, a 
measure for the number of 7r°s of which the decay photons were both detected is obtained by 
performing the integral over the thus obtained gaussian distribution. 

5 Results 

The method described above has been applied to the analysis of two sets of simulated data. 
Details concerning the simulation process will be given hereafter. 

The first set consisted of the data of 50,000 events, where the number of 7r°s produced per 
event was equal to 50. This dataset corresponds to peripheral events as detected in the WA98 
experiment. Of the 100 decay photons per event, about 16 were detected within the acceptance 
of the WA98 photon spectrometer [1]. 

In the fitting procedure of the T> distribution, only the bins corresponding to the range of 
90 to 180 MeV for the invariant mass were used. In fig. 4 the V distribution for this set of 
data is shown, together with the best fit V pre d- 
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Fig. 4. The T> distribution for the case of 50 7r°s/event. The full circles indicate the best fit T> pre( i. 

In fig. 5 the actual 7r° peak, as entered into the simulation, is presented as a dotted line, 
together with the reconstructed gaussian for the 7r° peak determined from the best fit for 
the reconstructed V distribution. The integral of the input 7r° peak over the fitted interval is 
57,479, whereas the result obtained with our reconstruction method amounts to 57,977. 
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Fig. 5. The input ir° peak (dotted line) and the reconstructed n° peak for the case of 50 7r°s/event. 

The second dataset consisted of 100,000 events. The number of 7r°s per event was 100, 
corresponding to medium central events as detected in WA98. Of the 200 decay photons per 
event, on average about 32 were detected. Fig. 6 shows the V distribution for this set of data, 
together with the best fit V pre d- 
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Fig. 6. The V distribution for the case of 100 7r°s/event. The full circles indicate the best fit V pre d- 

Fig. 7 shows the input 7r° peak (dotted line), together with the reconstructed gaussian for 
the 7T° peak. The integral of the input 7r° peak over the fitting interval is 229,712, whereas 
our reconstruction yields 224,023. 

It is seen that in both cases our reconstruction results are in very good agreement with the 
number of 7r°s as entered through the computer simulation. 
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Fig. 7. The input tt° peak (dotted line) and the reconstructed 7r° peak for the case of 100 7r°s/event. 
6 The computer simulation procedure 

In the computer simulation that is used to test the method, a fixed number of ir°s is created 
for each event. All of the 7r°s were decayed into a pair of photons. 

In order to obtain an exponential-like spectrum as observed experimentally for various parti- 
cle species [2,3,4], the ir° energies were distributed according to a Bose-Einstein distribution 
corresponding to a temperature of 200 MeV. 

The amount of 7r°s was chosen to match realistic conditions for high-energy heavy ion exper- 
iments. The number of accepted photons is comparable to the number of detected photons in 
peripheral or medium central events in the WA98 experiment. 

To mimic realistic experimental conditions, the directions of the photon momenta were 
compared to the acceptance angles of the WA98 lead-glass detector and only the photons 
within the acceptance were used in our analysis. In order to obtain a dataset that models the 
actual experimental data in a realistic way, the imperfections of the detector were simulated 
as well. The energy resolution of the detector, given by ^ = 0.02 + ^7= where all energies are 
in GeV, was used to obtain a smearing of the photon energies. Furthermore, we randomly 
discard 10% of all photons to simulate the effect of efficiency loss which in WA98 is due to 
charged particle vetoing. 



7 Summary 

The method introduced in this report allows accurate determination of the number of 7r°s 
via their two-photon decay channel in a range of 2a around the 7r° mass. The procedure 
does not invoke event-mixing techniques but is applied on an event-by-event basis using the 
difference between the total invariant-mass distribution (T) and a newly introduced " position 
swapped" distribution S. 

The method has been tested by means of computer simulations modelling peripheral and 
medium central events in Pb+Pb collisions at CERN-SPS energies. For these peripheral and 
medium central event samples it was seen that the ir° yield could be extracted with an ac- 
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curacy of 0.9% and 2.5%, respectively. As a final remark we would like to mention that the 
correlations of the T and S distributions have not been taken into account in the computation 
of the uncertainties of the "difference" distributions T>. However, it was seen that in the final 
results the effect of these correlations was negligible. 

In addition, we believe that this method, with some modifications, could be applied to the 
analysis of other two-particle decays, e.g. charm measurement via the D° — ► Ktt decay chan- 
nel. 
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